A cortical-inspired geometry for contour perception 
and motion integration 



In this paper we develop a geometrical model of functional ar- 
chitecture for the processing of spatio-temporal visual stimuli. The 
model arises from the properties of the receptive field linear dynamics 
of orientation and speed-selective cells in the visual cortex, that can 
be embedded in the definition of a geometry where the connectivity 
between points is driven by the contact structure of a 5D manifold. 
Then, we compute the stochastic kernels that are the approximations 
of two Fokker Planck operators associated to the geometry, and imple- 
ment them as facilitation patterns within a neural population activity 
model, in order to reproduce some psychophysiological findings about 
the perception of contours in motion and trajectories of points found 
in the literature. 

keywords: Visual cortex Lie groups Contact geometry Galilean group 
Cognitive neuroscience Spatio-temporal models 

1 Introduction 

The modular structure of the mammalian visual cortex has been discov- 
ered in the seventhies by the pioneeristic work of Hubel and Wiesel |21[|22| . 
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Abstract 




Many families of cells contribute to sampling and coding the stimulus im- 
age. Each family is sensible to a specific feature of the image: position, 
orientation, scale, color, curvature, velocity, stereo. Cells are spatially or- 
ganized in such a way that for every point (x,y) of the retinal plane there 
is an entire set of cells, each one sensitive to a particular instance of the 
considered feature, giving rise to the so-called hypercolumnar organization. 
Hypercolumnar organization and neural connectivity between hypercolumns 
constitute the functional architecture of the visual cortex, that is the deep 
machine underlying the processing of visual stimulus. 

To mathematically model the functional architecture of the visual cor- 
tex W.C. Hoffmann in [19| introduced a differential geometry framework in 
which the hypercolumnar organization was identified with a manifold earring 
a fiber bundle structure. Subsequent studies on the same line mainly focused 
on the set of simple cells, responsible for the detection of position and orien- 
tation. From a slightly different point of view Steven Zucker in |55| outlined 
that differential geometry instruments can be used to formalize the idea of 
boundary completion. Jean Petitot and Yannique Tondut described the set 
of simple cells as a contact structure |34], and introduced for the first time 
the notion of neurogeometry. Citti and Sarti proposed to represent the func- 
tional architecture in terms of Lie groups structures, which better describe 
the symmetry of the cortex [fjj. In their model, the hypercolumnar structure 
is described in terms of the Lie symmetries of the Euclidean group, equipped 
with the suitable sub-Riemannian metric. The integral curves of its gener- 
ating vector fields can also be considered as a mathematical representation 
of the association fields of Field, Hayes and Hess |14| . In this model contour 
completion is justified as a propagation in the sub-Riemannian setting |6j. 
In this structure they proved the relation between neural mechanisms and 
image completion (see also [41]). Analytical properties of their model was 
further studied by R. Hladky and Pauls [33]. Finally we recall the works 
of Duits, van Almsick, Franken, ter Haar Romeny |10| |11| |12| who pro- 
posed new models in different Lie groups, with many applications to image 
processing. 

In this paper we propose a mathematical model of cortical functional 
architecture for the processing of spatio-temporal visual information, that is 
compatible with both phenomenological experiments and neurophysiological 
findings. The adopted theoretical framework is the natural extension of the 
Lie symmetries approach proposed in [6,41 



Our starting point in Section [2] is a process of detection resulting from 
linear filtering with three dimensional Gabor functions with two spatial and 
one temporal dimensions, which have been proposed as a model of spatio- 



2 



temporal receptive profiles of primary visual cortex simple cells (see [9] , [7] ) . 
Spatio-temporal Gabor filters extend simple cells Gabor behavior as spatial 
filters 8p7 , that proved its usefulness for the classical task of edge detection, 



while its role for motion detection was already pointed out in [35]. The 
Gabor transform takes in input a moving image f(x,t), where x E M 2 are 
spatial coordinates and t E M. is time and provides in output a representation 
of the signal in the phase space 

f(x,t) -> F(q,s;p,u), 

with q £ R 2 and s E R representing spatio-temporal position and (p, v) 
representing, respectively, spatial and temporal frequency. 

Since we are mainly interested in spatio-temporal dynamical aspects, we 
will assign to temporal frequency the meaning of velocity v on the (x, t) 
coordinates. We will also select the subset of all detected features corre- 
sponding to a fixed value of \p\ . This will end up to be a 5D manifold, with 
a contact structure, induced by a normalization of the Liuoville form. 

In Section [3] we will outline that this constraint carries a notion of ad- 
missible curves [16] in a deterministic and a stochastic setting, allowing to 
compute the kernels connecting filters in the 5D manifold. The stochastic 
processes are completely determined by the described structures of the tan- 
gent space of the 5D manifold and they turn out to be described by the 
fundamental solution of a Fokker Planck equation 30 1 . In fact the processes 
contain diffusion in the fiber variables and transport along the remaining 
generators of admissible tangent directions, in the spirit of [29[ 52] . We will 
compute these kernels in two limit cases: the motion of a contour at a fixed 
time instant, and the motion of a point moving in time. The first one re- 
duces to the geometry of contours, with a notion of instantaneous velocity, 
the second one corresponds to point trajectories in time, and can be related 
to the subset of the Galilei group [43] on the plane. 

In Section [4] we will discuss the compatibility of the previously calculated 
kernels with psychophysiological findings reported in the recent literature. 
Then we will insert the connectivity kernels computed in Section [3] in a 
neural population activity model [13] , regarding them as cortical facilitation 
patterns. 

In section [5] we will use this population activity model equipped with the 
suitable connectivity kernels in two numerical simulations, comparing the 
results to recent phenomenological findings of the perception of contours in 
motion [36], and to fMRI measurements of cortical neural activity related 
to motion perception [54]. 
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Figure 1: A spatio-temporal receptive profile taken from the data set stud- 
ied in |9] and |7j (cell n.15 in the latter reference). In the top row, the 
sequence relative to the raw data, reconstructed using reverse correlation 
from electrophysiological recordings. In the bottom row, the 3-dimensional 
Gabor fit, visualized as 2D slices. 



2 The geometry of spatio-temporal dynamics 

In this section we extend an approach introduced in [6] that amounts to 
model each VI simple cell in terms of its receptive profile, to interpret its 
action as a Gabor filtering, and to introduce a geometry of the space com- 
patible with the properties of the output. In this paper we will consider 
each cell as sensitive to a local orientation and apparent velocity, that is the 
velocity orthogonal to a moving stimulus. The collected data lies on a five 
dimensional manifold Ai of space, time, orientation and velocity, a single 
cell being represented as a point on A4. The geometric structure of this 
manifold will be described in terms of a contact structure, that provides a 
constraint on the tangent space and on admissible connectivity among cells. 



2.1 Spatio-temporal receptive profiles 

It is known that the visual cortex decomposes the visual stimulus by mea- 
suring its local features. Local orientation and direction of movement have 
been the first visual features of neurons in VI that have been studied. Re- 
ceptive profiles (RPs) are descriptors of the linear filtering behavior of a cell 
and they can be reconstructed by processing electrophysiological record- 
ings [37] . It has been shown that the spatial characteristics of these RPs 
can be modeled by 2-dimensional Gabor functions [8 23 . 
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However a large class of cells shows a very specific space-time behavior in 
which the spatial phase of the RP changes gradually as a function of time [9] . 
Although many models have been proposed to reproduce these dynamics 
[T|9j , in [7] it has been shown that a 3-dimensional sum-of- inseparable Gabor 
model can fit very well experimental data of both separable and inseparable 
RPs (Fig. [TJ. Following this approach, we choose Gaussian Gabor filters 
centered at position q = (qi, q 2 ) G M 2 on the image plane, activated around 
time s G E, with spatial frequency p = (pi,P2) £ ^- 2 , temporal frequency 
v G R, spatial width a x G M + (circular gaussians) and temporal width 
a t G M+ 

_ \x-q\ 2 _ (t-s) 2 

^(x,t) = e i( - p - {x - q) - u{t - s)) e (1) 

where we have used the abbreviations z = (qi,q2,s) + i{j>\-,p%-,v) G C 3 and 
a = (a x , at) G M + x M + . Also note that the functions correspond to the 
propagation of a two dimensional plane wave within the activating window 
with phase velocity 

-H- 

The variable z = (qi, q2, s) + i(pi,P2, v) is canonically associated [15] to the 
phase space 

= {(qi,q2,s,pi,p2,f)} (3) 

endowed with the symplectic structure compatible with the complex struc- 
ture of C 3 of variable z, that is 

£1 = d\ = dpi A dqi + dp2 A dq2 — dv Ads 

where we have denoted with A the Liouville form 

A = pidqi + P2dq 2 - vds. (4) 

Since we are interested in analyzing the behavior of the receptive profiles 
at fixed spatial frequency, we introduce the reduced 1-form: 

uj = — = cosOdqi + sm.6dq2 — vds (5) 
\P\ 

where p = \p\ (cos 6, sin 9). This form is defined on the spatio-temporal phase 
space with fixed frequency, that is the 5-dimensional manifold 

M = R 2 x R+ x S 1 x R+ = {77 = ( qi ,q 2 ,s,e,v)}, 
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Figure 2: Visualization of the one- form oj with respect to a cell RP. In the 
left and center subplots are shown the isosurfaces of the RP already seen 
in Fig. [TJ In the left plot we visualize the reconstructed data, while in 
the middle plot we visualize its Gabor fit. The orientation of the yellow 
line within the fitted data represents the orientation of the wavefront of the 
Gabor filter and coincides with the one-form u = cos 6dq% + sin6dq2 — vds, 
or the dual vector X. In the right plot the semi-opaque red plane represents 
the horizontal tangent plane in the B? spatio-temporal space, generated by 
X$ and X\ (blue line). 
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and it is associated to every Gabor filter as shown in Fig. [2] At this level, 
time s is introduced as a base variable with the same role of q, with its 
dual v which has a similar role of 9, both being the engrafted variables with 
respect to time and space. To every point of M. corresponds univocally a 
Gabor filter whose parameters are the coordinates of the point itself. 

2.2 Admissible tangent space as constraint on the connec- 
tivity on Ai 

We will model the connectivity between points in the space M in terms of 
admissible tangent directions of M itself. From the geometric point of view, 
the presence of the 1-form ^ is equivalent to the choice of a vector field 
with the same coefficients as u, called dual vector (see Fig. [2]): 

X = (cosfl, sinfl, -u,0,0). 

The kernel of ui, denoted kerw (or cj = 0) is the space of vectors orthogonal 
to X. A basis of this space is constituted of the so-called horizontal or 
admissible vectors 

Xi = (- sin 9, cos 9, 0, 0, 0) , X 2 = (0, 0, 0, 1, 0) (Q) 
X 4 = (0,0,0,0,1) , X 5 = (v cos 9,v sin 0,1,0,0) 

so that 

kei uj = span{Xi, X2, X4, X^}, (7) 

that defines the horizontal tangent space. 

The only manifolds (curves or surfaces) admissible in this space are the 
ones whose tangent vectors are linear combinations of the horizontal vectors 
([6]). The neural connectivity between the receptive fields at different points 
in the M space will be defined in terms of these vectors in Section [3} 

For reader convenience we compute here explicitly all the non-zero com- 



7 



mutation relations between the vector fields El 

[Xi,X 2 ] = X 3 = (cos 6, sin 6, 0,0,0) (8) 
[^2,^3] = Xi 

[X 4 ,X 5 ] = X 3 (9) 
[X 2 ,X 5 ] = vXx. 

In Fig. |3] we have depicted the structure of the tangent fields Q. For 
visualization purposes, we show the fields in two different representations. 
In Fig. |3^l we visualize the structure restricted to the spatial and engrafted 
variables (q, 6, v), spanned by X\ and X 2 . The tilting of the planes is due to 
the non commutative relation ([8]). In Fig [3)3 we show the spatio-temporal 
structure restricted to the variables (23, s, v), where x% = q\ cos# + q 2 sin 0. 
Also in this figure the tilting of the planes spanned by X4 and X5 is due 
to the non vanishing commutator condition ([9]). It is worth noting that in 
this setting, the propagation is permitted only along the horizontal planes, 
while it is forbidden in their orthogonal direction X. Let us note that 
X3 = [Xi,X2] is linearly independent of the horizontal tangent space at 
every point. Hence (X\, X 2 , X4, X$) together with their commutators span 
the whole space at every point. This is the so called Hormander condition, 
which will have a crucial role in studying the property of the space. The 
same condition can be expressed in terms of properties of the form u, which 
is called contact form [16] since uj A dco is never zero being the volume form 
of the space. 

We also note that we can define on M a smooth composition law 
(q,s,e,v)Q(q',s',6',v') 

= (R e (q' + Qs') + q,s' + s, 6' + 9,v' + v) [W) 

where Rg is a counterclockwise rotation of an angle 6, that is such that the 
vector fields are left invariant with respect to this law. This implies that also 



To every vector field X = (ai, «2, a-a, 0,4, as) we can associate a directional derivative 
X = a 1 d qi + a 2 d qi + a 3 d s + a 4 <9g + a 5 d v 
with the same coefficients. Then we will call commutator of X and Y 

[X, Y]=XY - YX. 

We say that X and Y commute if [X, Y] = 0. Note that partial derivatives always 
commute, while directional derivatives in general do not. It is important to note that, 
even though the commutator is expressed formally as a second derivative, it is indeed 
a first derivative, so that it has an associated vector. Hence we can as well define the 
commutator between the vectors X and Y as the vector associated to XY — YX. 
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(a) (b) 

Figure 3: Visualization of the admissible tangent space of Ai, generated by 
the vectors X\, X2, X4, X5. In a) we show the projection in the (q, 9) space. 
Red planes are the sections of the horizontal planes, generated by X\ and 
Xi- Red arrows indicate the X\ direction, black arrows the vX% one. In 
b) we depict the projection in the (x$,s, v) space. The blue planes are the 
sections of the horizontal planes, spanned by X4 and X5 The black arrows 
represents the vXj, as before. 



the admissible curves of the structrue and its kernels will be invariant. M 
with this composition law can also be identified with a subset of the Galilei 
group (see Appendix). 



2.3 The output of the receptive profiles 

The output of the Gabor filters selects a set of points in M. corresponding to 
specific features of the image. Since the Gabor filters are always connected 
in terms of the vectors ([6]), also their output will inherit this structure, and 
will be concentrated around an admissible surface. 

The energy output of a cell with receptive profile in presence of a 
spatio-temporal stimulus f(x, t) is given by 



F(q,s,9,v) = \(r z ,f)\ 2 

dx dt e — CC 12 - - — <?)i cose+{x~q) 2 sing— u(t-s)) 



K3 



_ \x-q\ 2 _ (t-s) 2 

e °t f(x,t) 



(11) 



The analogous of ths formula (11) with purely spatial Gabor filters was 



studied by [6], where the filtering output F is a function of the variables 
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q and 9 alone. In that case, the output F to a stimulus f(x) takes its 
maximum around a value 9* = 6*{q) 

maxF(q,9) = F{q,9*) 
6 

where (cos 9*, sin 9*) identifies the orientation orthogonal to the level lines 
of f(x). As a consequence, in [6] it is proved that the level lines of / are 
lifted to curves admissible in the sense that their tangent vector lies in the 
kernel of the form 

W3 = cos 9dq\ + sin 0dq2 
which is a contact form that can be obtained by restricting ^ to the space 

M). 

Here we will take into account time, then the output F will be a function 
defined on Ai. In perfect analogy with the lower dimensional case, and 
denoting with (g, s(q)) a level set of the stimulus f(x, t), the output F takes 
its maximum around the values 9* = 8*(q, s(q)) and v* = v*(q, s(q)) 

max F(q,s(q), 9, v) = F(q,s(q),8*,v*) 

such that the vector (cos 9* , sin 9* , v*) is orthogonal to the level set of /. The 
vector (cos 9*, sin 8*) is orthogonal to the spatial level line, and the scalar 
v* represents the apparent velocity in this direction. 
The surface 

Z={(q,s(q),e*,v*)} 

is the 5D lifting of the level set of /, and the orthogonality condition implies 
that it is admissible, in the sense that its tangent vectors lie in the horizontal 
space, kernel of the form oj (see Fig. [3]). 

3 Curves and kernels of connectivity 

The lifted points of the spatio-temporal stimulus are connected by admissible 
integral curves. However we will see that not all admissible curves can be 
considered lifted ones, and we will describe which ones have this property. 
These curves will represent the association fields in space-time, analogously 
to the association fields of Field Hayes and Hess [l4| in the pure spatial case, 
and are reminiscent of the classical Gestalt concept of good continuation. We 
will discuss in Section [4] the compatibility of this structure with the cortical 
functional architecture. A probabilistic version of the connectivity field will 
be provided on the basis of Fokker Planck equations first introduced by D. 
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Figure 4: A schematical representation of the level set of a function / (in 
this case, an expanding contour) lifted as an admissible surface £ on the 5D 
structure. Red and black arrows indicate, respectively, the local direction of 
the vector fields X5 and X3 over two points in space-time. 

Mumford in the spatial case |29|, and interpreted as model of connectivity 
in Lie groups in |40| . Here we will extend this stochastic approach to the 
space-time contact structure. 

3.1 Generators of lifted curves 

We have already seen that the output of RP filtering is concentrated around 
admissible submanifolds. Nevertheless not all admissible submanifolds can 
be lifting of images, since lifting are graphs of the functions 9* and v* . For 
example the plane generated by the vectors X 2 and X4 cannot be recovered 
by lifting. Hence we will study the possible linear combinations of vector 
fields {A^}, which can be tangent to lifted level lines. An w-admissible curve 
7 C A4 is characterized by 7 G kera>, i.e. 

7 = a\Xi + a 2 X 2 + a 4 X 4 + a 5 X 5 

with at not necessarily constant. Lifted curves depend on the vectors X\ or 
X5, which are tangent to the base space (q, s), and their linear combination. 
To simplify the problem, we will consider separately the two limit cases of 
contour motion detected at a fixed time (05 = 0), and motion of a point in 
time (ai = 0). The first one is described by integral curves of the vector X\ 
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and of the generators of the engrafted variables X 2 and X4: 

73 = Xi + kX 2 + cX 4 . (12) 

These curves he in the section of the contact structure depicted in Figure 
[3^,, where k represents Euclidean curvature and c is the rate of change of 
local velocity along the curve. On the other hand, the motion of a point in 
time can be described as 

7T = X 5 + wX 2 + aX A (13) 

that are suitable to describe spatio-temporal trajectories of points. The 
coefficient w of the direction X 2 is the angular velocity, and the coefficient 
a of X4 is the tangential acceleration. 



3.2 Curves and kernels for contours in motion 

3.2.1 Moving contours as deterministic integral curves 

We consider now the geometry of moving contours at a fixed time. These are 
generated by the integral curves (12) so to satisfy the system of equations 

f j(t) = X^ + kX^ + cX^) 
I 7(0) = £0 

with variable coefficients k, c. In particular, due to condition ^ the set of 
points reacheable by piecewise constant integral curves of this type starting 
from a fixed point £0 = (Qo, so,6q,vo) is the space 

M so ={(q,s ,e,v)} (15) 

where so is fixed, so the space M. so can be identified with the space of points 

£ = (?,M). 

We also denote with £o(£o) the set of points reached by the fan of curves 
solution to the system (14) with constant coefficients (k,c), depicted in Fig. 

El 



3.2.2 Stochastic kernel 



Let us consider trajectories ( 14 ) with the additional property that the coeffi- 



cients k and c are allowed to vary at any point in a equidistributed way. This 



corresponds to substitute (14) with the vector- valued stochastic process 



(h; 

7(0) 



Xi( 7 )dt + X 2 {~i)dW l + X 4 ( 7 )dW 2 
£0 



(16) 
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Figure 5: The subset of horizontal curves with constant coefficients pertain- 
ing to So(^o) by varying parameters k and c, with initial velocity set to 0. 
Left plot: projections over the variables (<?i, <Z2> with different values for 
k, and same non-null value for c. Right plot: projections over the variables 
(<7i ,Q2,v) relative to some negative and positive values of c. The curves with 
k = (blue ones in the right plot) are straight lines, as the direction of 
spatial propagation never changes. 



where W = (WijWjs) is a two dimensional Brownian motion. The distri- 
bution of these stochastic curves is mostly concentrated around the surface 
So (£(})• The density of points reached by this stochastic kernel is then a can- 
didate to implement the mechanism of association fields (see Fig. |6j. This 
approach generalizes the approach of random paths introduced in |29|[40|[52] 
for the SE{2) problem. 

If we call po(£ 5 i|£o) 0) the density of points of M. So reached at the value 



t of the arc parameter by the sample paths of the process (16), then po can 
be obtained as the fundamental solution 

(dt + L )po(^ ,0) = 5^-^)S(t) 

where Lq is the Fokker-Planck operator 

L = X 1 - k 2 X\ - a 2 X\ (17) 

containing a diffusion over the fiber variables and v and a drift in the base 
variables, and with X{ and Xf we denote first and second order derivatives in 
direction Xj. Since Lq contains a set of vector fields that generates the whole 
tangent space of Ai so , by Hormander theorem this operator is hypoelliptic 
[20] . This means that po is non-null in all Ai So for any t > 0, even if the 
operator Lq contains only 3 linearly independent fields. Moreover, since 
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c. ^ 1, ^^^^^ q. ' 5 

(a) (b) 

Figure 6: Horizontal curves and stochastic kernels for the integration of 
contours in motion. Left: the isosurface plot of the kernel To (isovalue: 
0.002max(To)) is superimposed over the subsets of £o(£o) shown in Fig. 
[5j Right: the kernel projections of To relative to the variables a) (qi,q2,0) 
and b) (qi,q%,v) are plotted in gray under the projections of the curves 
(in yellow). This kernel was obtained by calculating the evolution of 10 6 
stochastic paths. 




it defines the Fokker-Planck equation for the stochastic process (16), po is 
indeed a (conditional) probability density on A4 S0 , evolving with the arc 
parameter t. 

In order to characterize each point of -M so in terms of the density of 



paths (16) that reach it, independently of the value of the arc parameter, 
we need a notion corresponding to the fan Eo(£o)j expressed in terms of 
Pq. The density of points reached at any value of the arc parameter by the 



stochastic dynamics ( 16 ) is given by 

ro(£|£o)= I Po(C,t\Co,0)dt . (18) 

This derived quantity is actually the fundamental solution of the operator 
Lq, so explicitly we have 

Xx - { K 2 d 2 9 + a 2 d 2 v ) )r„(£|£o) = S(Z - Co)- (19) 

The fundamental solution ro(^|^o) can be directly implemented in terms 



of random paths beginning in £o by solving numerically the system (16). 
An example is shown in Fig. |6j where an isosurface plot of the kernel is 
plotted over the same integral curves of Fig. [5] From the figure it can be 
seen that the probability density is concentrated around the surface So (Co) 
and decays rapidly away from it. 
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3.3 Curves and kernels for point trajectories 

3.3.1 Point trajectories as deterministic integral curves 



The trajectory-type curves introduced in (13) are solutions 7 
the system of ordinary differential equations 

= X 5 ( 7 ) + wX 2 (-f) + aX A { 1 ) 
= Vo 




M to 



(20) 



for given initial point 770 = (^iq, (?2o> sq, 9o, vo). In general the coefficients 
w,a need not to be constant, but we can have a local approximation of 
any curve in a neighborhood of the starting point if we consider the case of 
constant coefficients. In this model case when w and a are not zero, (20) is 
explicitly solved by 



cos((9 - 4>)) + —tsmOit) 
w 

sin(6> - 4>)) - —tcos9(t) 
w 





= Qio + g(cos(9(t) 


«&(*) 


= Q20 + g(sm(9(t) 


< s(t) 


= s + t 


9(t) 


= 9 + wt 


> v(t) 


= vo + at 



where g 



— and d> = arctan , and a reasonable choice is to set 

so = 0, in order to synchronize the arc parameter t with the time parameter 
s. 

The fan £(770) of such curves is depicted in Fig. [7j Each of them de- 
scribes a motion on the plane q that, for small times, consists approximately 
of arcs of circles with radius g, while for sufficiently large times corresponds 
to enlarging spirals around a slightly moving center. Their instantaneous 



acceleration is given by 



gential acceleration^] along direction X3 



Rg\ , so that a corresponds to the tan- 

K vw, 



a = qi cos 9 + (J2 sin 9 
2 By direct computation, a is the time derivative of the modulus of the velocity 
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Figure 7: Projections over the variables (qi, q 2 , s) of the fan of curves with 
constant coefficients, for different values of parameters w and a and a non- 
null initial velocity. Different values of w are associated to different curve 
colors, while gray is used for projections over two dimensional planes. Left 
plot: curves with different values of w, and a fixed non-null value for a. 
Right plot: curves with different values of a. 



while their curvature is 

glgg ~ gggi _ w 
(qi 2 + q2 2 ) 3 * ~ v ' 

In particular, when w = we obtain straight lines along the Oq direction, 

while for a = we obtain circular trajectories of radius — . 

w 

3.3.2 Stochastic kernel 



Similarly to what we have done in section 3.2.2, let us consider the vector- 
valued stochastic process 

7(0) = m [ j 

where W = (Wi, W 2 ) is a two dimensional brownian motion. 

As before, we see that the density p(rj, t\r)o,0) of points of M. reached at 



the value t of the arc parameter by the sample paths of the process (21 ), is 
the fundamental solution of the equation 

(d t + L)p( V ,t\r]o,Q) = S(r}-r} )5(t) 
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(a) (b) 

Figure 8: Horizontal curves and stochastic kernels for the integration of 
motion, with a non-null initial velocity value. As before, different values of 
k are associated to different curve colors. Left: the isosurface plot of the 
kernel T (isovalue: O.OOlmax(r)) is superimposed over the curves showed 
in Fig. [7J Right: the projections of T relative to the variables (qi,q2,s) 
are plotted under the projections of the curves (in yellow). This kernel was 
obtained by calculating the evolution of 10 6 stochastic paths. 



where L is the Fokker-Planck operator 

L = -a 2 X\ - k 2 X\ + X 5 . (22) 

Integrating the density p with respect to the arc parameter, we derive the 
fundamental solution of the operator L. Explicitly we have 

d s - U 2 8 2 + a 2 8 2 v - vX 3 ) )T( V \ m ) = 5(rj - m ) (23) 



hence this is now a propagation in the physical time s. 

Similarly to what we have done in the previous section, we can obtain 
T(rj\rjo) at each point by solving numerically the system of stochastic differ- 



ential equations (21). An isosurface plot of the kernel is shown in Fig. |8j 
over the same integral curves of Fig. [7] and again we see the concentration 
around the fan £(770). 

The fundamental solution (23) will be concretely applied in Section [5] as 
a facilitation field, toward the task of motion integration of trajectories. 



4 Neural propagation of boundaries and trajecto- 



ries 



The curves (14) and (20) can be related to well-defined perceptual mech- 



anisms. More precisely, we propose to consider them as association fields, 
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in the sense of |14| , devoted to the integration of contour in motion and 
trajectories. The perceptual tasks of subjective boundary completion and 
motion integration have been widely studied by both psychologists and phys- 
iologists, and their underlying physiological explanation continues to be an 
open issue of discussion. The perceptual bias towards collinear stimuli has 
classically been associated to the long-range horizontal connections linking 
cells in VI sharing similar preferences in stimulus orientation. This special- 
ized form of intra-striate connectivity pattern is found across many species, 
including cats, tree shrews and macaques [5j[24j, the main difference be- 
ing the specificity and the spatial extent of the connections. Furthermore, 
axons seem to follow the retinotopic cortical map anisotropically, with the 
axis of anisotropy being related to the orientation tuning of the originating 
cell |4j. These connections have already been modeled in [6] by means of a 
sub-Riemannian diffusion process over the orientation space R 2 x S 1 '. 

Similarly to what happens for the integration of merely spatial visual 
information, the brain is also capable to easily predict stimulus trajecto- 
ries |44|, and to group together boundary elements sharing similar motion 



or apparent motion paths 25 , 36 . One possible explanation for these ef- 
fects could be the existence of specialized facilitatory networks linking cells 
anisotropically and coherently with their preference in velocity and axis of 
motion direction. Supporting these assumptions, it has been found that the 
neural preferences in direction of movement are also structurally mapped 
in the cortical surface, with nearby neurons being tuned for similar motion 



direction 50 , and it has been shown that excitatory horizontal connections 
in the VI of the ferret are strictly iso-direction tuned |38|. Furthermore, it 
is known that also extra-striate area MT/V5 is retinotopically organized, 
its horizontal connectivity pattern being highly structured, with connec- 
tions reaching columns of cells tuned for similar orientation and direction 
preference anisotropically and asymmetrically [2 28 . Moreover, striate and 
extra-striate cortical areas seem to cooperate, and surround modulation in 
VI can be given by the connectivity patterns implemented in both areas by 
means of fast feedforward and feedback inter- areal projections |3j. 

Regarding the physiological correlates of motion integration, it was shown 
that the motion of an occluded object trajectory is significantly represented 
in the human brain by the same visual areas that process real motion |31| . 
Moreover, a recent study based on electrophysiological recordings of the VI 
of tree shrews showed some non-linear neural behaviors that are coherent 
with the phenomenological dynamics of motion integration |54|. It is in- 
deed a general assumption that some cortical area in the visual cortex is 
responsible for predicting future motion, a possible implementation being 
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a specialized connectivity for spatio-temporal trajectory facilitation, that is 
different from the one in VI responsible for contour integration 18||45[|4"7"||51| . 



Here, we do not speculate on the exact physiological origin of these psy- 
chophysiological findings. However, we want to show that the connectivity 
kernels Tq and T arising from the geometry defined in Section [2] is capable to 
reproduce qualitatively some of the effects reported on the works that have 
been previously cited. We will embed this connectivity in the neural popu- 
lation activity model described in the next paragraph. Then, in Section [5] 
we will simulate the response of cortical visual neurons to artificial stimuli, 
comparing the results with some psychophysiological findings reported in 
the recent literature. 

4.1 Modeling neural activity 

The state of a population of cells can be characterized by a real-valued 
activity variable, which depends on the interaction of the feedforward input 
at different points, due to the cortical connectivity. The first population 
activity models are due to Wilson and Cowan [53], and Ermentrout and 
Cowan 1 131: 



da(rj, t) 



-a(v, t) + S(c f J rfa|C)a(C, t)d( + F(r?)) (24) 



where a is the neural activity of the population, F is the feedforward input 



(11), r is the facilitation kernel, ct is the facilitation strength and S is the 
sigmoidal function 

Sir) = • (25) 



In the stationary case a first order approximation of the solution of (24) is 



F F ( V ) = s(c f f T( V \ V )S(F( V ))d( + F(rjfj , (26) 

that is the activity formula that we will use in the experiments of Section 
[5] Let's explicitly note that the term 

F T ( V ) = S(F( V )) (27) 

represents the mean neural extra-cellular activity in response to a stimulus. 

The geometry of the functional architectures is contained in the kernel 
r, taking into account the deep structures of the connectivity space. Then 
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the overall probability of activation can be obtained by convolution of the 
activity with the kernel T, and the term 



P(rj) = J r(r ? |r ? )5(F T (7 ? ))dr ? (28) 

is the facilitation pattern resulting from the contribution of horizontal cortico- 
cortical or feedback inter-areal connectivities. 

In the model of Ermentrout-Cowan only position and orientation were 
considered and symmetry properties were imposed to the facilitation pat- 
terns. In [40j, for the same features it was proposed to choose as a kernel 
r the fundamental solution of a Fokker Planck equation deduced from the 
Euclidean symmetries. In the next section, two numerical simulations will 



be performed using equation (26). In the first one, we will model the prop- 



agation of boundaries in motion using the connectivity kernel (18). In the 



second experiment we will model instead the propagation of trajectories, 



using the 5D kernel (23) 



5 Numerical simulations 

In this section some numerical simulations will be performed to test the reli- 
ability of the kernels computed in Section [3] when introduced in the activity 



model (26) and the results will be compared to some psychophysical and 



physiological experiments. 

5.1 The feedforward and extracellular activity in response 
to a stimulus 

For the subsequent numerical simulations, we measured in every discrete 
point Zj = (xj,yj,tj) of a stimulus / of size n x x n y x n t , its local energy of 
orientation and speed by convolving the input image sequence with a pre- 
determined bank of Gabor filters of fixed spatial frequency centered at the 



points rn = (xi,yi,ti,9i,Vi), thus discretizing Eq. (11) so to have 

N 

Ffa) = \y~^ e - i (pi( x j- x i)+P2(yj-yi)-M t j- t i)) 

3=1 (29) 

_ \^j-^i\ 2 + \Vj-Vi\ 2 _ (tj-'i) 2 2 



where p\ and p2 are the component frequencies \p\ cos 6i and \p\ sin^j. 
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frame 1 frame 14 frame 26 frame 39 frame 51 frame 64 




Figure 9: The stimulus used in Experiment 1. 



We chose the spatio-temporal frequencies and the maximum velocity 
value v m represented in the Gabor filter set depending on the stimulus, so 
that given the couple (|p|,w m ) the filters have a maximum temporal fre- 
quency of u m = \p\v m . We normalized the filters so that the response could 
range from (in regions with no changes in luminance) to 1 (square plane 
waves sharing the filter parameters \p\ and fj), corresponding to a normal- 
ization of the stimulus contrast. The Gabor scale parameters will follow the 
relations 

2.57T 

(30) 



Ox — 4 |p| 



2u„ 



whose meaning is to approximately have 2.5 spatial subregions under the 
Gabor Gaussian envelope, and a variable Vi that allows the filter with maxi- 
mum velocity to cover one wavelength within the Gabor active time interval 



(see I7U9 37 for the physiological justifications) 



The thresholded feedforward input Ft is then computed following Eq. 
(25), representing the mean neural extra-cellular activity in response to the 
stimulus. In the subsequent analyses, the values to which we set the param- 
eters \i and j3 in the above stated formulas will be explicitly specified. It is 
worth noting that these parameters are not due to physiological data, due 
to the difficulty in finding some reference in the literature, but are chosen 
as to have a computationally reasonable number of non-null measurements 
F T . 



5.2 Experiment 1 - Contours in motion 

For the first simulation we use a dashed circle in linear motion moving up- 
wards within the visual space. The image size is 200 x 200 pixels and the 
image sequence is composed of 64 frames. The circle segments are approx- 
imately 2 pixels wide and the circle is moving with a uniform speed of 0.5 
pixels per frame (see Fig. [9]) . To get a measurement of local spatio-temporal 
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features, we convolve the stimulus with a set of 3-dimensional Gabor func- 



tions with (11). The spatial frequency parameter is set to \p\ = n/2, so that 
the Gabor moving subregions match the width of the segments. The max- 
imum local velocity represented in the filter set is v m = 1 pixel per frame, 
while spatial and temporal scale parameters are calculated accordingly with 



(30). After having convolved the stimulus with the Gabors, we model the 



neural extra-cellular activity by using (25) with parameters fi = 10, f3 = 0.5, 
obtaining Ft- 

We select the 4-dimensional subset corresponding to the 32nd frame 
FT v (£,i) = {FT(rji)\ti = 32}, thus neglecting behaviors over time, where 
£i = (xi,yi,6i,Vi). This is because the model given by the stochastic model 
of connectivity distribution Tq is stationary, even if it is relying on spatio- 
temporal information, and is compatible with the dimensionality of Ai SQ in 
(fl5). 



Following the assumptions that the measurements Ft v (£) can model the 
output of a cortical direction- and speed-selective cell (or of a neural popula- 
tion that are selective to the same visual features), the overall continuation 
probability can be obtained with a discretized version of the cortico-cortical 



facilitation pattern (28) with the kernel Tq: 

^ = E r °&l^(£), ( 31 ) 

where Fo(£i|£^) is calculated for every fiber vector on the same dis- 

cretized domain M. SQ of the lifted stimulus, using the stochastic approach 
described in Section [3.3.21 

It is worth noting that the parameters k and a governing the diffusion 
process are related to the maximum spatio-temporal curvature of an illu- 
sory contour. Here, we have chosen the values for the diffusion coefficients 
manually in order to have the best results, because of the lack of reference 
physiological data. In particular, the parameters were set so that at the final 
value taken by the arc parameter of the stochastic paths, the mean square 
displacements of the fiber variables (6,v) are (ir,v m /2). 

Finally, the population activity is computed by 

F F (Z i ) = S(F&) + c f P(£ i )), (32) 



that is a discretized version of (26) where cj is a coefficient governing the 



total strength of the excitatory connections. In Fig. 10 are shown the re- 
sults of the simulation. In the top row we have plotted a iso-level surface 
of the normalized and thresholded output measurements Ft v . The process 
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Figure 10: Results of Experiment 1: in the left subplots, measurements of 
local direction of motion 9; in the right subplots, measurements of local 
velocity v. The theoretical values for a circle translating at a speed of 
0.5 pixels per frame are shown with a red line, a) Isosurface plot of Ft 
(isovalue: 0.2). The information about local orientation and velocity is 
correctly retrieved, but activity remains clustered within disjoint regions 
of the domain, b) The modeled extra-cellular activity Fp after horizontal 
propagation on the submanifold M SQ (facilitation strength: cj = 40). 



correctly lifts the stimulus around its theoretical values in the domain Ai so . 
In the bottom row, the iso-level surface of the global cortical activity Fp is 
shown. Note that continuation is performed within the whole manifold M SQ , 
propagating the local speed and direction cues at the end of the segments, 
and interpolating the data to define a 4-dimensional set that carries infor- 
mation even in the visual space between the segments, where no changes 
in contrast were originally present in the stimulus. These results are co- 
herent with recent psychophysiological findings, where it is stated that a 
global shape in motion is better perceived if local velocity changes smoothly 
along its contour |36j, and that the integration of the motion of a partially 
occluded object is facilitated when its visible contours define closed config- 
urations p7| . 



5.3 Experiment 2 - Motion integration 

The stimuli that we will use in this simulation will be several instances of an 
object in motion along a certain direction, that disappears at a given time 
position t\, and reappears at t2 with a direction of motion changed by A9 in 
a position that is coherent to a piecewise continuous trajectory (Fig. 11 ). It 
is a well documented fact that humans tend to perceive the two trajectories 
as pertaining to a single unit just for small values of AT = ti — t\ and A8, 
an effect that is commonly referred to as motion integration. In particular, 
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frame 1 frame 7 frame 13 frame 20 frame 26 frame 32 




time 

Figure 11: An example stimulus sequence, and the paradigm used to gener- 
ate the stimuli in Experiment 2. 



we know from many psychophysiological experiments that the chances of 
detecting a straight or curvilinear trajectory in noise increase with stimulus 
duration and with direction coherence [46j[49] . 

To do so, we created multiple instances from the stimulus paradigm 



described above, and showed in Fig. 11 assigning different values to the 



parameters AT and A9. It is worth noting that we are using an elongated 
object whose axis of motion is always orthogonal to its major axis of ec- 
centricity. The choice of this particular stimulus is naturally explained by 
the geometry, as the points £ Ai that are connected by its continuation 
property implicitly define a local orientation and a direction of motion that 
is orthogonal to that orientation. The Gabor filtering of a different kind of 
stimulus (for example a moving dot) would indeed measure high responses 
also for those velocity projections that are coherent with the real axis of 
motion, due to the motion streak effect [17] , requiring a different neural 
model to detect unambiguously the direction of motion. Thus, even if the 
majority of the psychophysiological experiments use moving dots as stimuli, 
the validity of the following experiment is not influenced, as exploring the 
exact physiological implementation of motion integration is out of the scope 
of this paper. 

The image sequences to process are 51 x 51 pixels wide and are composed 
of 102 frames. The value of eccentricity of the moving ellipsoidal object is set 
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(a) (b) 

Figure 12: A comparison between the spatio-temporal continuation capa- 
bilities of the model, applied to stimulus instances with different values of 
AO. The plot shows three isosurfaces of maxw „) Fp (isovalues 0.9, 0.5, 0.1 
in red, yellow an green): a) AT = 12, AO = f , b) AT = 12, AO = ff . The 
probabilistic continuation given by the geometry in the configuration with 
the smallest change in direction of motion allows a smooth trajectory com- 
pletion. The same is not true for the stimulus with an higher value of A9: 
even if a weak response connecting and interpolating the two parts of the 
stimulus is still present, the strongest facilitation component do not deviate 
from the early straight path of the object. 



to 2, and its velocity is set to 0.5 pixels per frame. Each stimulus instance is 
uniquely identified by the couple (AT, AO). We convolve each sequence with 
a set of Gabor filters with v m = 1 and p set to match width of the object's 
minor axis, and again, we threshold the output according to (25) using 
parameters [i = 20, (3 = 0.7, obtaining the set of measurements Ft- In the 
previous simulation we have taken a temporal slice of this output, and then 
we have propagated the activity using the connectivity kernel IV Similarly, 
here we will propagate the activity present in Ft, but without discarding 
time, thus convolving it with the 5-dimensional kernel V defined on the same 
discrete domain M of the lifted stimulus, using a formula similar to (31). 
The corresponding facilitation pattern P(rji) is then obtained and the total 
population activity Tp(rjj) is computed following (32) with Cf = 20. 

The results obtained by processing two stimulus configurations are showed 



in Fig. 12, where only the central portion of the full temporal domain is 



plotted, to highlight the different effects that activity propagation has in the 
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Figure 13: A direct comparison between the electrophysiological recordings 
made in [54] and the results obtained with our model for the same stimulus 
parameters: in both cases AT = 0, while A8 is equal to | (a) and | 
(b). To compare responses, we integrate over x and y and depict three 
isosurfaces with isovalues taken as 0.9, 0.5, 0.1 (in red, yellow and green 
respectively) times the maximum value among the outputs showed. Note 
that, disregarding the weaker physiological responses to opposite directions, 
the non-linear facilitating behavior is very well reproduced by our model. 
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two cases. When the change in direction of motion of the object is under 
a certain threshold, the trajectory is completed smoothly, yielding a strong 
activation in the temporal interval when the object is not present in the 
scene. If the value of AO increases, however, the spatio-temporal curvature 
of the optimal subjective trajectory progressively becomes too high to be 
completed by the connectivity, even if a weak response linking the two parts 
of the stimulus could yet be present. 

It is worth noting that due to the thresholding stage, this subjective 
interpolation effect is strongly non-linear, as it cannot be explained by the 
sum of the responses to the first (i, < t\) and second (t{ > tz) part of 



the stimulus alone (see Fig. 13). This is coherent with the findings in 
the work of Wu et al, where electrophysiological experiments recorded a 
similar non-linear behavior of trajectory interpolation for small values of 
AO [54] . In that work, the cause of this effect is left unexplained. Even 
if the stimulus paradigms are slightly different (they use a field of moving 
dots with a null value of AT) , the result of our simulations allow us to make 
some speculations. 

A good qualitative description of the effects that the modeled excitatory 



connectivity has on stimulus response can be found in Fig. 14 where for 
every stimulus configuration we plot the difference 

F f ac = F F (S 3 )-F F (S 1 )-F F (S 2 ) (33) 

where S3 is the full stimulus and Si and S 2 represent the first (51 (ti > 
t\) = 0) and the second (S2(ti < t 2 ) = 0) part of the stimulus. The visu- 
alization of the output Ff ac highlights the role that a trajectory-specialized 
cortical connectivity could have in performing tasks of motion integration. 
The spatio-temporal extension of Ff ac over Ai gets smaller for higher values 
of AT and AO. The decay of the facilitation effect is coherent with the 
observations made in [47] (where AO = 0) and [54] (where AT = 0), even 
if the experiments are methodologically different. Moreover, even if, coher- 
ently with our results, some psychophysiological experiments showed that 



a broken trajectory in noise is easily detectable 42,48 , as far as we know 
little has been done to explore the effect of changing the duration AT of the 
temporal occluding gap. 

5.4 Discussion 

Regarding the implementation of our model, we would like to highlight func- 
tional meaning of the parameter couple (n,a), driving the diffusion along 
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the fiber variables (0, v) when calculating the kernels Tq and T. In particu- 
lar, the parameter k seems to be strongly related to the maximum perceived 
curvature of illusory contours and trajectories, while a sets the maximum 
rate of change of local velocity along admissible subjective contours (tra- 
jectories). It would be interesting to try to fine-tune the parameters of 
the model in order to reproduce quantitatively as precisely as possible the 
psychophysiological findings found in literature. We aim to carry out this 
analysis in a future paper. 

It is worth noting that the physiological correlates of the first simulation 
are not well documented and, as far as we know, no one has yet studied 
the neural activity in cortical regions responding to subjective contours in 
motion. We suppose that it would be a difficult issue to address, as multiple 
cortical layers, as VI, V4 and MT/V5, may be involved [26]. Recently, some 
significative results have been obtained with static illusory contours, using 
promising electrophysiological techniques [32] . Those kinds of experiments, 
targeted to the detection of illusory motion contours, could give additional 
clues about the neural computation that governs the influence of spatio- 
temporal features in the detection of moving shapes and boundaries. 

As discussed in the beginning of Section |4j a possible physiological im- 
plementation of the geometry used in the second numerical simulation, could 
be a trajectory-specialized cortico-cortical connectivity between neurons in 
the higher visual areas. The advantages of having such an excitatory con- 
nectivity implemented in the visual cortex could be very important in per- 
forming complex cognitive tasks. For example, we assume that the brain 
could use the improvements in contour detection with respect to the only 
SE{2) related functional architecture, due to an additional information, that 
of velocity, that indeed in practical situation is a feature that is coherent 
on objects. The influence that the connectivity arising from our geometrical 
model has in the processing of visual tasks such as spatio-temporal grouping 
or segmentation will be the subject of a future paper. 

6 Conclusions 

In this paper we proposed a model of cortical functional architecture for 
the processing of spatio-temporal visual information. Motion features are 
detected first by simple and complex cells with RP modelled by 2D+time 
Gabor filters. Linear filtering with such a profiles liftes the visual stimu- 
lus from the base space R 3 to the phase space M 6 comprising spatial and 
temporal frequencies, in which a Liouville form can be defined. 
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Figure 14: The non-linear effect of the facilitation by varying the stimulus 
parameters AT and AO. Each subplot shows the isosurface of Ff ac relative 



to different stimulus instances (isolevels calculated as in Fig. 13). The non- 



linear effect of the facilitating kernel V remains appreciable just for small 
values of AT and AO, rapidly decaying as the value of the characteristic 
stimulus parameters increases. 



We defined a 5D phase space with fixed frequency by taking a reduc- 
tion of the previous differential form. Then, exploiting the commutation 
properties of the horizontal basis, we regarded the tangent space of the 
contact manifold A4 as a possible constraint acting on the connectivity 
between points, giving the definition of admissible integral, or horizontal, 
curves. Possible linear combinations of the horizontal basis that are com- 
patible with the definition of admissible curve have been studied, in order 
to model the possible lifting of the visual stimuli as association fields, in the 



sense of 14 . We considered the corresponding deterministic integral curves 
for two modeling limit cases: contours in motion and trajectories of a point 
in motion. 

Then we considered horizontal stochastic paths, i.e. trajectories of points 
that always move along the tangent space of A4 and are allowed to change 
the value associated to the fiber variables (0, v) in a random, equidistributed 
way. We have seen that the resulting integrations over the arc parameter of 
the paths coincide with the kernels of the Fokker Planck operators defined 
on the geometries. Both deterministic curves and stochastic kernels inherit 
the symmetries of the non associative reduced Galilean group as described 
in Appendix. 

After a discussion about the compatibility of our theoretical frame- 
work with various phenomenological and psychophysiological findings on 
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visual perception and cognition reported in the literature, we introduced 
the stochastic kernels as facilitation inducers in a neural population activity 
model. 

Suitable numerical simulations have been carried out, by processing pre- 
determined artificial stimuli with the neural population activity model pre- 
viously described, showing that the Fokker Planck stochastic kernels endow 
the model with the capability of completion and continuation of contours 
in motion and trajectories of points, coherently with the phenomenological 
experiments of Rainville [36] and Wu et al |54| . 

In conclusion, results have shown that the proposed functional geometry 
is compatible with existing psychophysical and physiological experiments, 
even if a complete knowledge about the effective neural implementation 
needs supplementary empirical data. 



Appendix: Relations with the Galileian group 

The geometry of (M,cj) can be related with the classical Galilei group of 
motions on the plane, when choosing an appropriate local frame of reference. 

The two dimensional Galilei group is the semidirect product Lie group 
G = (Rjj x R s ) x (Rl x Si) defined by the composition law [43] 

(<z, s ,M)-(«V,0V) (u) 

= {R e q' + us' + q,s' + s, 9' + 9, R e u' + u) { ' 

where (q, s) provides space-time translations, Re is a counterclockwise pla- 
nar rotation of an angle 9 and u = (^) is the velocity vector responsible of 
the purely Galilean transformation, or boost, expressed in Cartesian coor- 
dinates. 

The manifold Ai can be identified with a subset of this group by iden- 
tifying the direction of velocity with the angle of the rotation. Setting the 
velocity in polar coordinates 




= v cos (f 
= v sin ip 



then we can identify A4 with the set of points 

M = {(q,s,9,Rgv)} 

that is, using the identification 9 = (p. Since the angle 9 is interpreted as 
the local direction of a boundary, then this amounts to select the apparent 
velocity orthogonal to it. 
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This manifold is not invariant under the Galilei composition, but it is 
invariant with respect to the smooth composition law ( 10 ). This composition 
law is not associative, but possesses a trivial neutral element e = and each 
rj € M possesses a left inverse element, i.e. t]~ 1l £ Ai such that rj~ lL Qrj = 

(q, s, d, v)- lL = (-R- e (q - Qs), -s, -6, -v) 

and for these reasons we can consider it a reduced Galilei non associative 
group p9] . 

In particular, standard invariance arguments apply to the solution to 
( 23 ) , that possesses the symmetry 

r(7/|jfo) = r(?? " lL ©#) Vi),i)oeM. 



Analogously, the fan 2(770) of integral curves (20) with constant coefficients 
starting from a point r/o can be obtained from the one starting from by 

S(lo) = 1o lL 0S(O) . 
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